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ABSTRACT 

Cosmic shear constrains cosmology by exploiting the apparent alignments of pairs of galaxies due 
to gravitational lensing by intervening mass clumps. However galaxies may become (intrinsically) 
aligned with each other, and with nearby mass clumps, during their formation. This effect needs to 
be disentangled from the cosmic shear signal to place constraints on cosmology. We use the linear 
intrinsic alignment model as a base and compare it to an alternative model and data. If intrinsic 
alignments are ignored then the dark energy equation of state is biased by ^ 50 per cent. We examine 
how the number of tomographic redshift bins affects uncertainties on cosmological parameters and 
find that when intrinsic alignments are included two or more times as many bins are required to obtain 
80 per cent of the available information. We investigate how the degradation in the dark energy figure 
of merit depends on the photometric redshift scatter. Previous studies have shown that lensing does 
not place stringent requirements on the photometric redshift uncertainty, so long as the uncertainty 
is well known. However, if intrinsic alignments are included the requirements become a factor of 
three tighter. These results are quite insensitive to the fraction of catastrophic outliers, assuming 
that this fraction is well known. We show the effect of uncertainties in photometric redshift bias and 
scatter. Finally we quantify how priors on the intrinsic alignment model would improve dark energy 
constraints. 

Subject headings: large-scale structure of Universe - cosmological parameters - surveys 
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1. INTRODUCTION 

The tidal gravitational field of density inhomogeneities 
in the universe distorts the images of distant galax- 
ies. This so-called 'cosmic shear' results in correlations 
in the observed ellipticities of the distant galaxies, a 
signal which depends upon the geom etry of the uni- 
verse and the matter power spectrurn (iBlandford et aD 
fl99lt iMiralda-Escudg fl99l iKaised [199^) . The first 
detection s of cosmi c shea r made in 2000 ([Bacon et"al] 
2000: K aiser et all [2OOOI: IVan Waerbeke et al.l l2000r 
Wittman et al.ll2000[ ) demonstrated its value as a cosmo- 



logical tool. Future generations of multi-color imaging 
surveys will cover thousands of square degrees and have 
the potential to measure the dark matter power spectrum 
with unprecedented precision in three dimensions at low 
redshift which is not possible with the CMB. This is im- 
portant because dark energy only starts to dominate at 
low redshift. Ultimately, cosmic shear data can improve 
constraints on the dark energy equati on of state from the 
CMB by over an order of mag nitudellxJ (l2002f ). 

The great promise and exactitude of cosmic shear has 
necessitated the design of instruments expressly geared 
towards measurement of the tiny lensing-induced dis- 
tortions. It has also motivated improvements in tech- 
niques to account for changes in galaxy shapes due to 
the atmosphere and telescope optics. Furthermore it has 
prompted careful consideration of any potential cosmo- 
logical contaminants of the cosmic shear signal. Intrinsic 
alignments of galaxies are a potential contaminant and 
fall into two categories. The first is intrinsic-intrinsic 
galaxy alignments (II correlations) , which may arise dur- 
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ing the galaxy formation process since n eighboring galax- 
ies res ide in a similar tidal field (e.g. ICrittenden et al.l 
(|2001h ). The second, related, effect is a cross-term be- 
tween intrinsic ellipticity and c osmic shear (GI correla- 
tions, iHirata and Seliala ()2004D ). whereby the intrinsic 
shape of a galaxy is correlated with the surrounding den- 
sity field, which in turn contributes to the lensing distor- 
tion of more distant galaxies. The net effect of this is 
an induced anti-correlation between galaxy ellipticities, 
leading to a suppression of the total measured signal. 

Intrinsic galaxy alignment has been subject to nu- 
mero us analytic, numer i cal a n d observationa l stud ies 
[e.g. [C roft and Mctzlcr' (2000'); 'Heavens et al.' (2000 
ICrittend e n et al, (20 01); Mack ev ct al. (2002); Jin;. 
II2002D: iHevmans et all (^2004a^: 'Hirata et all (12001 



IBrid le and Abdallal (120071 )]. iBrown et al.l (|2003[) and 
Hev mans et al. I (I2004a|[ detected an II s i gnal in the Su- 
perCOSMOS data. iMandelbaum et all (|2006l ) used in 
excess of a quarter of a million spectroscopic galax- 
ies from SDSS to obtain constraints on intrinsic align- 
ments, with no detection of an II signal. Since II 
correlations only operate when galaxies are physically 
close, this contaminant is relatively straightforward to 
deal with by down-weightin g or remov ing those pairs 
if redshifts are well known (iKing and Sc hneider 200^; 
IHevmans and Heavens! 120031: iTakada and White .2004) 
or by fitting parameterized models relying on the dif- 
ferent behavior of lensin g and II in redshift space 
(jKing and Schneidedl20M l. 

The first observational detection of a large-scale 
density-galaxy ellipticity correlation was made by 
IMandelbaum et al.l (|2006( l. using the same SDSS sam- 
ple noted above; their detection is robust and is present 
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on scales up to 60 Mpc. They estimate that the am- 
phtude of the GI correlation could cause existing deep 
surveys to underestimate the linear amplitude of fluctu- 
ations by as much as ~ 20%. The GI signal is dominated 
by the brightest galaxies, possibly due to these being 
BCGs (bri ghtest cluster galaxie s) aligned with the cluster 
ellipticitv. iHirata et al.l |2007f ) perform a more detailed 
characterization of this effect, including a higher redshift 
sample of LRGs (luminous red galaxies). They estimate 
that results on erg from future cosmic shear surveys may 
be bi ased by around 10 per c ent. Using N-body simula- 
tions, iHevmans et al.l ()2006af l estimate that in a survey 
with median depth ~ 1, the GI signal can contribute 
up to 10% of t he lensing signal on scales up to 20 arcmin. 

As noted bv lHirata and SelialJ (|2004D . unlike II corre- 
lations, GI correlations are not localized in redshift, and 
hence their removal is more complex. Given that the 
GI term has been observati o nally detected we cannot 
ignore it. iHirata and Selialj ()2004l ) suggested that the 
dependence of the GI signal on redshift could be used 
to distinguish it from cosmic shear. King ()2005f ) showed 
that projecting the signal into a set of template functions 
would enable the lensing, II and GI signals to be isolated, 
again harnessing their different dependence on redshift. 

In this paper we aim to remove the intrinsic alignment 
contamination of cosmic shear by simultaneously mea- 
suring cosmological and intrinsic alignment parameters 
from shear power spectra. The increased number of fit- 
ted parameters may degrade constraints on cosmologi- 
cal parameters. Here we focus on the requirements this 
places on photometric redshift quality. 

The planned future imaging surveys will observe at 
least hundreds of millions of galaxies. It will not be 
possible to obtain spectroscopic redshifts for them all. 
Therefore it is planned to rely on 'photometric redshifts': 
estimates calculated from galaxy luminosities in several 
observing bands. The resulting redshift accuracy is criti- 
cally dependent on the number of observing bands used. 
Equally important for measurement of cosmological pa- 
rameters is the existence of a representative sample with 
spectroscopic redshifts which allows the redshift accu- 
racy to be quantified. It is important to plan now for 
these future observations since decisions are being made 
about which observing bands to use. Furthermore it may 
be necessary to mount a dedicated observing campaign 
if more spectroscopic redshifts are needed than are cur- 
rently planned. 

In practice additional information may be included in 
the fit to help offset the degradation due to using many 
free parameters to encapsulate the intrinsic alignment 
mo del. This may come f rom m easu rements suc h as those 
by iMandelbau m et all ()2006[ ) and IHirata et al.l (|2007l ). 
Any additional information on at least the functional 
form of the intrinsic alignment contributions from the- 
ory would also be a great help. This information could 
be included by applying priors to the intrinsic alignment 
power spectra in the simultaneous analysis. We inves- 
tigate the size of prior required to improve constraints 
relative to the self-calibration regime. 

In section [2] we introduce our fiducial intrinsic align- 
ment model. Section [3] shows how dark energy con- 
straints are affected by varying parameters within this 
model. In section |4] we show the requirements this places 
on photometric redshift quality. Finally in section [5] we 



consider the effect of priors on photometric redshift and 
intrinsic alignment parameters. 

2. INTRINSIC ALIGNMENT MODEL 

In this paper we attempt to simultaneously measure 
the intrinsic alignment and cosmological model from cos- 
mic shear power spectra. We therefore have to parame- 
terize the sources of the intrinsic alignments in some suf- 
ficiently flexible yet physically reasonable way. We also 
need to define a fiducial model for intrinsic alignments 
to put into our simulation of the future observations. 

We assume a catalogue of source galaxy positions 
and shears which is divided up into a number of bins 
in redshift, with a number per unit comoving distance 
ni{x) for bin i. For convenience this is normalized as 
J ni{x)dx = 1. The lensing efficiency function for lens- 
ing a mass at redshift for source redshift bin i may 
then be written as 

^ c Jq Xs 

where Xd a-nd Xs Sive the comoving distances from the 
observer to the deflector and source respectively. We 
assume a flat universe throughout. Qm is the matter 
density parameter, Hq is the Hubble constant and c is 
the speed of fight. 

As detailed in IHirata and SeliaS (|2004f ) the shear 
power spectra between redshift bins i and j come from 
the sum of three terms 

where the first term is the usual gravitational lensing 
contribution, the second term arises from the intrinsic 
alignments of physically close galaxies and the final term 
arises from the intrinsic alignment of a galaxy with a 
mass which lenses a more distant galaxy. 

In this paper we consider only E mode power spectra, 
since the lensing and intrinsic alignment contributions to 
the B mode power spectra are both very small (for the 
model considered). 

These contributions can be written in terms of under- 
lying power spectra as 

c.l) = r^^%^P.(fc;x)rfx (3) 

Jo X 

/•°° nj(x)wj(x) p X, fr^ 

'^('^'>Jf^ ^7nK;x)ax (5) 

where k = i/x- The first power spectrum Ps{k;x) is 
simply the (non-linear) matter power spectrum at the 
defl ector redshift. For s implic ity in this paper we use 
the P eacock and Doddsl (|1996f ) method to modify a lMal 
(1999) power spectrum (containing baryon supression 
but no wiggles). 

The remaining two power spectra are not well known. 
We use perturbations around the linear alignment model 
normalized approximately to data. We detail our as- 
sumptions and compare them to data below. 

2.1. Fiducial intrinsic-intrinsic term (II) 
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Fig. 1. 



Dotted line: The best fit to the SDSS data using 



the HRH* model. Dot-dashed line: The prediction from the linear 
alignment model roughly normalized to IHirata et al.l (|2004| ) , who 
normalize to SuperCOSMOS. Solid line: The prediction from the 
linear alignment model if a non-linear matter power spectrum is 
used in place of the linear theory matter power spectrum. The solid 
line is closer to the H RH* model and we use th is model through- 
out. The error bars of lMandelbaurn" encompass all the 
models. 

Tlte intrinsic-intrinsic term describes tfie extent to 
wliicli nearby galaxies are aligned with each other, for 
example due to common tidal forces during their for- 
mation. For this paper we simplify t h e line ar align- 
ment mod el presented i n ICatelan et al.l ()2001[ ) and de- 
ve loped in lHirata etal.l (|2004j) to use only the first term 
of lHirata~et al . (20Q3) Equation 16, written 
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(6) 



where Ci is a normalization constant, and Pg^^{k) is 
the linear theory matter power spectrum. D{z) = 
(1 + z)D{z) where D{z) is the growth factor normal- 
ized to unity at the present day. p is the mean mat- 
ter density of the Universe as a function of redshift. 
We estimate the yalue o f Ci by matc hing to the power 
spectra in Figure 2 of iHirata et al.l (p0b4il and esti- 
mate Ci = 5 X 10~^^{h^MQ/Mpc-^)-^. IHirata et all 
(12004) chose their Ci by comparison with SuperCOS- 
MOS lBrowneTall (I200I . 

The term we ignore is roughly an order of magnitude 
s maller than the ter m we include. This can be seen 
in lHirata et al.l (|2004j) by examining the B-mode H power 
spectrum which is of the same order as the term we ig- 
nore. Indeed from the figures it can be seen that the 
B-mode H power spectrum is typically an order of mag- 
nitude below the E-mode H power spectrum. 

The linear alignment model might be a reasonable ap- 
proximation on large scales, but clearly makes a number 
of simplifying assumptions on sm all scales. Inspired by 
the remark in Hirat a et al.l ([2007) (Section 8) we attempt 
to make a more realistic model by inserting the full non- 
linear matter power spectrum Ps^k) intolH] 



(7) 



We refer to this throughout as the non-linear l inear 
alignment model. As noted in IHirata et al.l ()2007f ) this 



"model" has no grounding in theory. It might not even 
be closer to the truth: galaxies may become more or less 
aligned due to non-linear interactions. However we find 
this "model" matches slightly better to other models and 
data and therefore we use it as our default throughout. 

The HRH* model for the i ntrinsic alignraent cor rela- 
tion function was proposed bylHeavens et al.l (|2000b[ ) and 
improved by iHeymans et al.l ( 20041)1 ). This model pre- 



dicts the average intrinsic alignment on the sky between 
two galaxies, as a function o f the three dimensional sepa- 
ration between the galaxie s. IHeymans et al . (2004b) and 
iMandelbaum et aTl (|2006( l placed constraints on the am- 
plitude of this model using COMBO-17 and SDSS data 
respectively. 

These correlation functions may be projected along the 
line of sight to produce quantities, w++ and Wxx, which 
are closely related to the II power spectrum in Eq. O by 



P;yi{k) = 27r J w++{rp)Jo{krp) Vpdrp 



(8) 



where rp is the separation between pairs of galaxies in 
the plane of the sky, and throughout Jn{F) are Bessel 
functions of the first kind. Here Wj^.^ is the correlation 
of galaxy shapes using the component of shear aligned 
along the line joining the two galaxies. There is another 
correlation function Wxx which comes from the shear 
components at 45 degrees to the line joining the two 
galaxies. In the above we have followed HRH* and as- 
sumed that these two correlation functions a re th e same 
(u)++ = Wxx) and thus simplified IHirata et al.l ()2004l ) 
Equation 10 to write Eq. [H 
The HRH* model prediction for w++ is given by 



w++irp) 
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1 + 1- 



-7gs 



1 



1 + {r/B)- 



■ drii 



... 

where TZ is the factor to convert measured ellipticity to 
shear, r|| is the separation between galaxy pairs along the 

line of sight, and r — -|- . The galaxy clustering 

correlation function is described by tq and 7^^, and A 

a nd B are free param et ers of the HRH* model. 

IMandelbaum et "aTl (|2006( ) follow IHeymans et al.l 
(|2004b[ ) by fitting A to data and using jgg = 1.8, 
B = l/i-^Mpc, and tq 5.25 /i"^ Mpc. For their 
method TZ ~ 0.87. We calculate w++ using these pa- 
rameters and the be st fit to SDSS L3-L6 of A = 1.8. We 
follow Section 4.3 of lMandelbaurn et al.l (|2006( ) by using 
a range of integration of < r\\/{h~^Mpc) < 60h~^ 
Mpc. The result is shown by the dotted line in Fig. [TJ 
Our X axis range covers the range of data used for the 
fit. 

We use the linear alignment model to predict by 
inverting Eq. [51 



w++irp) = 



1 

2^ 



P^i {k)jQ{krp) kdk. 



(10) 



We calculate this at a redshift of z — 0.12 to match 
the mean redshift of the SDSS sample. This is shown 
in Fig. [1] by the dot-dashed line, which uses the linear 
theory matter power spectrum (Eq. [6] in Eq. llOp . The 
solid line shows the result when the non-linear theory 
matter power spectrum is used (Eq. [7]in Eq. llOp . 
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Fig. 2. — The dot-dashed line shows the linear alignment model. 
The solid line shows our fiducial model, which uses the non-linear 
matter power spectrum in the linear alignmen t model. Data point s 
taken from the analysis of SDSS in Figure 1 of 'Hirata et al.l II2007I '). 
Circles are for blue galaxies and crosses for red galaxies. Note 
that only rp > 11.9h~^ Mpc (rp > 4.7/i~^ Mpc) were used in 
their most (least) conservative fits (grey shading). Upper L4, the 
largest sample. Lower L6, the brightest sample. The upper two 
lines show the linear alignment model predictions multiplied by a 
factor of ten. 



It can be seen that the Hnear ahgnment model predic- 
tions have a similar amplitude to the HRH* fit. This is 
not surprising since the amplitude of t h e line ar alignment 
model was taken from 'Hirata et al.| (|2004f) who made 
a fit to the iHev mans ct al. (2004 b) d ata which is not 
far from the iMandelbaum et al.l (j2006D result. The lin- 
ear alignment model is lower than the HRH* model on 
small scales and larger on scales 3/i~^Mpc < < 6Qh^^ 
Mpc. However the HRH* model was developed for small 
scale correlations therefore the discrepancy with the lin- 
ear alignment model on large scales is not a concern (it 
was based on an approximation to the density correla- 
tion function that only applies on small scales). Both the 
linear alignment model and the HRH * model fit comfort - 
ably within the error bars of Mandc lbaum et ahl (|2006D . 
We note that the non-linear linear alignment model is 
closer to the HRH* model. 



2.2. Fiducial shear-intrinsic term(GI) 



As shown in IHirata et al.l ()2004D , the linear alignment 
model predicts for the GI underlying power spectrum 
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\k). (11) 

where Ci is the same number as in Eq. [6l We also define 
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where Pg is the full non-linear matter power spec- 
trum which may help to encapsulate non-linear ef- 
fe cts on smaller scal e s. F or example it was shown 
in iBridle and Abdallal ()2007| ) that even if there were no 
tidal interactions between galaxies, there will still be a 
shear-intrinsic correlation on small scales if galaxies are 
aligned with their own dark matter halos. 

We compare this to the SDSS results by IHirata et al.l 
((2007I) . who measure the real space galaxy-shape correla- 
tion function projected along the line of sight Wg+ which 
is defined by 



C3+(rp,H)dH 



(13) 



where ^g+(rp,H) = (7+}(rp.n) is the radial shear of one 
galaxy relative to another, averaged over all pairs of 
galaxies of a se paration r„ on the s ky and H along the 
line of sight (see IHirata et alJ ()2007l ) equation 10 for the 
estim ator they use) . 

The IHirata et al.l (|2007f ) results for luminosity bins L4 
and L6 are shown in Fig. |2l We show their results for L4 
because it has the largest number of galaxies and there- 
fore will be most typical of future cosmic shear surveys. 
We show results for L6 (which contains the most lumi- 
nous galaxies) because there is a detection of the signal in 
this band and we would like to compare the shape of the 
linear alignment model prediction with the observations. 

We can more easily predict theoretically the mass- 
shape correlation function wg-^- which can be related to 
the corresponding power spectrum by 

ws+{rp) = J Ps,=yi{k)J2{krp) kdk (14) 



(|Hirata et all l2007l Equation 23). IHirata et all (|2007D 
convert from ws+ to Wg+ by multiplying by the galaxy- 
mass bias of the sample which we take from their Table 
5. We use the effective redshift va lues from thi s table 
and the same cosmological model as IHirata et al.l (|2007l ) 
to calculate Ps,^' and convert it to Wg+ shown in Fig. |2l 

The L4 data points fit very nicely over the fiducial 
linear alignment model predictions. Note that the dark 
(light) gre y shaded reg ions block over the data points 
that lHirijita et al l (j2007t ) did not use for their most (least) 
conservative cuts when fitting models. It is not trivial 
that these data agree with the model. The amplitude 
C'l in the linear alignment model was chosen to match 
observations of the II power spectrum. The agreement 
with the GI data is some support for the linear alignment 
model which links the two. 

The L6 data points are much above our fiducial linear 
alignment model predictions. This is because we are aim- 
ing to use parameters correct for the majority of the pop- 
ulation, whereas L6 is a very bright sample. We multiply 
the predictions by a factor of 10 to compare their shapes 
with the data. We see that there is a slightly better fit 
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Fig. 3. — A selection of shear cross power spectra for our fiducial survey divided into ten tomographic redshift bins. We assume a 
photometric redshift scatter of &z = 0.05(1 + z) and zero catastrophic outliers. Solid line: lensing shear (GG) term. Dashed line: intrinsic 
shear (II) term. Dotted line: shear-intrinsic (GI) cross term (absolute values shown). Dot-dashed line: total. Galaxies are divided into 
tomographic redshift bins with equal numbers of galaxies in each bin. The median redshifts of the tomographic redshift bins 1 to 10 are 
(0.30, 0.49, 0.62, 0.73, 0.84, 0.96, 1.1, 1.2, 1.4, 1.9) respectively. 

D UNEt and LSST^ the 'sha llow' survey described 
in lAmara and RefregieH (|2006D . This covers 20,000 
square degrees to a depth of Zm = 0.9 with 35 galaxies 
per square arcmin ute which can be used for sh ear mea- 
surements. As in lAmara and Refregieil (l20p6f) we use 
the redshift distribution given bv lSmail et all (119941 ) us- 
ing a — 2 and f3 ~ 1.5. By default we divide this survey 
up into twenty redshift bins, binned to have an equal 
number of source galaxies in each bin. As discussed be- 
low we use dz to parameterize the photometric redshift 
error, and by default use Sz = 0.05. 

To show the size of the intrinsic alignment contribu- 
tions we plot in Fig. [3] some of the power spectra for a 
survey divided into 10 source redshift bins. The GI and 
II contributions to the lensing power spectra are typically 
about an order of magnitude lower than the GG contri- 
bution. For illustration we show only the results for a 
high redshift subset of the bins. However, for the auto 
correlation of the lowest redshift bin the II signal is the 
largest contribution to the total power spectrum. When 
cross correlating the lowest redshift bin with the 2nd and 
3rd lowest bins, the GI term can become larger than the 
GG term and the total power spectrum negative. 

Despite the small size of the intrinsic alignment contri- 
bution, we must note that the effect on the power spec- 
tra of the cosmological parameters, including the dark 
energy equation of state, can be very smal l and therefore 
these contributions are very important. iHuterer et al.l 

t The Dark UNi verse Explorer (DUNE) 
|http : //w ww, duiie-misslon.net/ 1 

? The Large Synopticl3k}^Survey (LSST) http : //www ■ Isst . org] 



when using the non-linear matter power spectrum in the 
linear alignment model equation (Eq. [T^ . It was effec- 
tively this comparison which originally led iHirata et al.l 
( |2007f ) to suggest using the non-linear power spectrum. 

On scales less than 4.7/i~l Mpc the L4 results appear 
to agree slightly better with the linear theory predictions 
than those from the non-linear matter power spectrum. 
However the error bars are large and these data points 
are not used in the model fitting by Hirata et al because 
the bias is non-linear on these scales. We expect the non- 
linear alignment model to be more difficult to disentangle 
from the cosmic shear signal because it has a similar 
shape, as a function of scale. We show in Fig. [7] the 
results for both models and find them to be quite similar 
in their photometric redshift requirements. 



2.3. Implications for our fiducial survey 

For the remainder of this paper we use as o ur fidu- 
cial c osmological model the best fit to WMAP3 (ISpergell 
l200(]f) . This has a normalized Hubble constant h = 
i?o/100kms~^Mpc~^ — 0.73, total matter density — 
0.23, baryon density ili,h^ = 0.022, fluctuation ampli- 
tude normalization erg = 0.74 and primordial fluctuation 
spectral index ris = 0.95. We use fiducial dark energy 
parameters wq = —0.95 and Wa — where the dark 
energy equation of state as a fu nction of scale factor is 
given by w(a) = w q + (1 — a)wa (jChevallier and Polarskil 
l200lt iLinded [20031) . We use a maximum angular power 
spectrum multipole Zmax = 2 x lO'' throughout. 

Throughout we use a fiducial survey similar to 
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TABLE 1 

Percentage biases on cosmological parameters 
when wrongly ignoring intrinsic alignments. 
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Fig. 4. — Illustration of the default spatial flexibility we allow in the intrinsic alignment power spectra. Lines are as in the previous 
figure but the intrinsic alignment power spectra have been multiplied by a function with five bins in k which are interpolated linearly in 
the log as detailed in the text. We have arbitrarily set the bin values to (n— i) ~ ^' *° show the freedom. 

()2006t ) give an approximate expression (equation 24 of 
their paper) for the dependence of the lensing power spec- 
tra on cosmological parameters which includes C^*^ oc 
trl'^lwol"'^"'^. Therefore an increase in the fluctuation 
normalization parameter erg by 1 per cent increases the 
power spectra by around 3 per cent. A change in w of 1 
per cent (e.g. w = —1 to w — —0.99) typically changes 
the lensing power spectrum by about 0.3 per cent (de- 
crease). With future surveys we are aiming to measure 
the dark energy equation of state to of order 1 per cent. 

We simulate power spectra that include intrinsic align- 
ment contributions but find the best fit cosmological pa- 
rameters that would be obtained if intrinsic alignments 
were mistakenly ignored (GG). The purpose is to assess 
the importance of bothering to consider intrinsic align- 
ments for future studies. We make the simplistic assump- 
tion that changes in the cosmic shear power spectra are 
linear in the cosmological parameters. If this assumption 
were correct then Fisher matrix uncertainties would also 
be correct, so it is a useful investigation to compare with 
predicted cosmological parameter uncertainties. 

Given the numbers above, it is not surprising that if we 
fit to the lensing power spectra mistakenly ignoring the 
GI and II contributions we mis-estimate the fluctuation 
amplitude parameter cg by ^ 5%, in agreement with 
previous studies (see the upper section of Table 1 for 
details). From the numbers above we expect the effect on 
the dark energy equation of state to be about a factor of 
10 larger. Indeed we find that the dark energy equation 
of state is mis-estimated by around 50 per cent for a 
non-tomographic (2D) analysis. 



This mis-estimate only applies if all cosmological pa- 



NOTE. — Upper section: for a non-tomographic (2D) 
analysis. Lower section: for 10 bin tomographic analysis. 
Columns show: the percentage change in erg when fixing 
all other cosmological parameters; the absolute change in 
WQ when fixing all other parameters; and for wo when 
marginalising over the other 6 cosmological parameters 
used in this paper. The results can be quite different 
when marginalising over other parameters due to move- 
ment along degeneracy directions to better fit the intrinsic 
alignment signal. 

rameters are known except the dark energy equation of 
state. In fact we find that if several cosmological param- 
eters are fitted simultaneously the mis-estimate becomes 
even worse, as near-degeneracy directions are exploited 
in an attempt to match the shape of the GG-I-GI-I-II 
power spectra using the GG spectra alone. Therefore if 
our fiducial model is roughly correct it would be a dis- 
aster to ignore intrinsic alignments in studies aiming to 
measure the dark energy equation of state. 

Table ?? shows that the final overestimate of the dark 
energy equation of state is 15% for a 2D non-tomographic 
analysis, and over 400% for a 10 bin tomography anal- 
ysis. Note that these numbers are only indicative, since 
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they assume that the hkeUhood surface is Gaussian with 
a constant curvature, which wiU not be a good approxi- 
mation so far from the fiducial modeL 

The simplest way we can hope to remove the II con- 
tribution is to use only cross-power spectra and to ig- 
nore the auto-correlations. Fig. [3] reminds us that more 
than just the auto-correlations may need to be removed if 
the photometric redshifts are not perfect and many red- 
shift bins are used. For 10 photometric redshift bins and 
6z = 0.05 the II contribution is non-negligible even when 
cross correlating non- neighboring bins (e.g. i = 6, j = 8). 
Further, due to the photometric redshift errors the GI 
contribution is non-negligible even for auto-correlations. 

2.4. Parameterizing intrinsic alignment uncertainties 

On obtaining shear power spectra from new data we 
will attempt to estimate cosmological parameters. How- 
ever, since the linear alignment model is approximate we 
cannot rely on it being correct. It seems unlikely that 
the galaxy formation process will one day be understood 
so well that we can use it to perfectly remove intrinsic 
alignments. We can hope to measure the contaminant 
from analyses suc h as th ose of lMandelbaum et al.l ()2006f ) 
and iHirata et al.l (j2007l ) . However it will still be nec- 
essary to use a parameterized model to fit to the data 
points, and there will remain some uncertainty on the 
model parameters due to the error bars on the data 
points. 

For this paper we parameterize perturbations around 
the linear alignment model and show the effect of fitting 
these free parameters simultaneously with the cosmo- 
logical parameters. This follows the app roach indicated 
in section 7 of the technical appendix of lAlbrecht et al.1 
|2006). Our most basic perturbation is to assume some 
unknown amplitude and redshift dependence, which is 
different for each of II and GI: 



X 



1 



1 



zo 



Px{k;x) (15) 



where X is 7^ or S, 7^ . 

However this assumes that the scale dependence of the 
models is well known, whereas we have seen that there 
is still great uncertainty (e.g. compare the HRH model 
with the linear alignment model in Fig.[T]). Further, the 
scale will likely deviate in a different way from the linear 
alignment model at different redshifts. In principle we 
have two unknown 2d functions Q^j (k; x) and Qs,j' {k; x) 
multiplying each of the power spectra for the II and GI 
terms: 

P'r{k-,X) = Qx{k;x)Px'''ik;x)- (16) 

To parameterize the Qx using a finite number of pa- 
rameters we use n bins in k and m bins in redshift and 
linearly interpolate in the logs of k, (l + z) and Qx- That 
is, 

HQxik, X))=KZ Bf^ +{l-K)Z Bf^^,^^ 
+K (1 - Z) Bf^^_,,) + {1-K){1- Z) S(^+i)(^.+i)(17) 
where 

[ln(A:) - ln(fc,)] 



K- 



ki < k < ki+i 



[\n{k,+i)~\n{h)] 

[ln(l + z)- ln(l + z,)] 
[ln(l + zj+i)-ln(l + z,)] 



(18) 



(19) 



and i runs from to n and j runs from to m. Be- 
cause we are interpolating on the log of a multiplicative 
function setting all the values to zero makes the mul- 
tiplicative function everywhere unity. 

We use different free parameters i?^ for II than for 
GI. However we use the same number of bins n, m for 
both of II and GI. We set fcg and A:„+i to be at the 
edges of our k integration ranges (fco = lO^^'/iMpc"'^, 
fc„ = 2 X lO^/iMpc"^) and fix i?oj = i?(„+i)j = for aU j 
i.e. the multiplicative function goes to unity at low and 
high k. ki to kn are spaced linearly in the log in a smaller 
range in k. We use ki = 0.1/iMpc~^ and kn = 2/iMpc~^ 
because this spans from near the start of the non-linear 
regime to a point which has less effect on our angular 
power spectra. 

The basic effect on the angular power spectra is illus- 
trated in Fig. m Here we have used no evolution in red- 
shift (7X = 0, m = 1) and set = (1, —1, 1, —1, 1) 
to highlight the degree of flexibility and the £ values cor- 
responding to the edges of the bins. 

By default we use one bin in redshift (m = 1) which 
means that we assume the same power law redshift evo- 
lution in each k bin. By default we use n = 5 bins in 
k. We explore the effect of these defaults below. For nu- 
merical stability we apply a wide prior on each B value 
of 10. We see in Section [5] that this is sufficiently wide 
that it would give the same answer as a wider prior. 

3. DARK ENERGY CONSTRAINTS 

We now consider the constraints on dark energy that 
can be placed after marginalising over other cosmological 
and intrinsic alignment parameters. We give our results 
in terms of the Figure of Merit (F oM) advocated by th e 
Dark Energy Task Force report (|Albrecht et al.ll200(l . 
This is the inverse of the area of the 95 per cent contour 
in the wq, Wa plane. We marginalise with flat priors 
over firm 1^81 and Us assuming a flat Universe. 
Note that the WMAP3 cosmological parameters give a 
smaller FoM compared to most other sets of parameters, 
since they have lower values for ag and n^. For 
acceptable values of these parameters the absolute value 
of the FoM can be a factor of two or more higher for 
the survey parameters we assume. However this roughly 
scales the FoM the same for all subsidiary parameter 
values we consider. 

We use a Fisher matrix appro ach to quantify un c ertain - 
ties. We use the equations in iTakada and Jainl (|2004l ). 
We add the GI and II power spectra to the GG power 
spectra and insert this as the simulated total power spec- 
tra. When calculating the Fisher matrix we need the co- 
variance matrix between all the power spectra. For the 
cosmic variance part we use the sum of all the power 
spectra. Because the GI term is negative this sum can in 
principle be zero. Therefore we checked the effect of us- 
ing only the lensing (GG) term in the covariance matrix. 
This gives almost the same results, except the figures 
of merit are slightly bigger, particularly for the GG-f-II 
results due to the decreased error bars on ignoring the 
positive II contribution. 

In principle the linear alignment model predictions de- 
pend on cosmology via the growth factor, matter density 
and matter power spectrum. In principle this depen- 
dence could be used to constrain all cosmological param- 
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Fig. 5. — Top row: Assuming perfect photometric redshifts [S^ = 0). Bottom row: More realistic photometric redshifts [5^ = 0.05). 
Left column: Including just the II terms. Middle column: Including just the Gl terms. Right column: Including both II and Gl terms. 
Solid lines: figure of merit as a function of number of kbins. Lines from top to bottom within one panel: the number of bins in the redshift 
direction increases 1, 2, 5. Dotted line: FoM for GG alone. Dashed line: using the linear alignment model with the linear theory matter 
power spectrum, instead of the non-linear theory matter power spectrum. In every case a free amplitude parameter and unknown power 
law evolution was marginalized over (allowed to be different for each of Gl and II). The default intrinsic alignment model we use for the 
remaining figures is shown by a circle and our maximally flexible model is marked by a cross. 

more bins in wavenumber (x-axis of panel) and redshift 
(upper to lower lines) both degrade the constraints by a 
similar but small amount. 

The lower left panel of Fig. \E\ shows that in practice 
even if P;yi(k; x) has few free parameters it will not help 
to constrain cosmology better than if intrinsic alignments 
were negligible. This is because with realistic photomet- 
ric redshifts it is not possible to adequately disentangle 
the GG and II contributions to the shear power spec- 
tra. However the amount of degradation is still rela- 
tively small and flat with increasing number of bins in 
wavenumber. 

The upper and lower middle panels of Fig. [5] show the 
result of marginalising over an unknown Gl power spec- 
trum in the absence of any II contribution. The degra- 
dation is more severe than for the case with II alone. 
We attribute this to the greater similarity between GG 
and Gl angular power spectra (e.g. Fig [3]). Therefore 
it is more difficult to disentangle the two contributions. 
Again, the degradation is most sensitive to the number of 
redshift bins m. For the most flexible model we consider 
the degradation is about a factor of two independent of 
the photometric redshift error. 

Including both II and Gl terms at the same time de- 
grades the dark energy constraints yet further. With 
zero photometric redshift errors the GG-I-II result is es- 
sentially the same as the GG alone result. We note that 
also the result for GG+GI+II is almost the same as that 
for GG-f Gl. Thus the II parameters are fitted adequately, 
presumably using the information from the autocorrela- 
tion power spectra. However when more realistic photo- 
metric redshifts are considered we see that the II term 
does cause a degradation. 
It is not physical to allow complete freedom in the 



eters if the intrinsic alignments alone could be observed. 
However, we do not believe sufficiently strongly in the lin- 
ear alignment model to use it in this way. We therefore 
calculate the intrinsic alignment power spectra P;yi (fc; x) 
and Ps,^i once, for our fiducial cosmological model, only. 
When we make Fisher steps in cosmological parameters 
we continue to use these intrinsic alignment power spec- 
tra for the fiducial model. This makes it impossible to 
constrain parameters such as as, and Hs through their 
effect on the linear alignment model. This seems sensi- 
ble to us, and in practice has little effect on our results. 
The intrinsic alignment angular power spectra then de- 
pend on cosmology only through the weightings in Eq.s|4] 
and [51 

First we consider the effect on the FoM of unknown in- 
trinsic alignment terms in the case of perfect photomet- 
ric redshifts {6z =0). If no II or Gl terms are included 
FoM = 38 for our fiducial survey. The upper left panel of 
Fig. m shows that when the II term is included we can ac- 
tually get an increase in FoM if sufficiently little freedom 
is given in the intrinsic alignment model! This is because 
the II angular power spectra (Eq. [5]) do depend on cos- 
mology even though we have fixed Pxji {k; x) to be inde- 
pendent of cosmology. Because we have here assumed 
perfect photometric redshifts then the autocorrelation 
power spectra measure very well the II angular power 
spectra separately from the GG angular power spectra. 
Features in the II angular power spectra subtend different 
angular scales on the sky and thus constrain cosmology 
through the standard ruler effect. The size of this effect 
will depend on the fiducial model used for example, if it 
is a power law with wavenumber then it would not help 
to constrain cosmology. This constraint weakens as more 
free parameters are used in the II power spectra. Adding 
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intrinsic alignment power spectra. There are vari- 
ous physical processes that cause the power spectra to 
vary with scale and redshift. For example, the pro- 
cesses t hat tidally align galaxie s with neigboring mass 
chimps (jHirata and Seli ak* 2004") may evolve differently 
to those which align galaxy light and mass, an d thus 
cause GI correlations (jBridle and Abdallal [2"007f) . It is 
very unlikely that there are so many physical processes 
that the power spectra could, for example, oscillate in 
the redshift direction 5 times between redshift 3 and 0. 
Similarly it is hard to think of physical mechanisms to 
cause ten changes in power law slope in the k direction. 
Our most flexible model has 52 free parameters for the GI 
underlying power spectrum and the same again for the II 
power spectrum. Considering this large number of free 
parameters it is perhaps surprising that the cosmologi- 
cal parameter constraints are not weakened still further. 
We take as our default model for the remainder of the 
paper n — 5, m = 1, which means we have freedom in 5 
bins in the wavenumber direction as illustrated in Fig. |3] 
and the same power law in the redshift direction for each 
wavenumber bin. We also show some results from our 
maximally free model with n = 10 and m = 5. 

We repeat all the calculations using the linear mat- 
ter power spectrum in the linear alignment model (Eq. [6] 
and Eq. [TT|) . shown by the dashed lines in Fig. [5l We 
see that most of our results are qualitatively similar, al- 
though the FoMs are generally lower. We attribute this 
to the difficulty of constraining the parameters of the 
intrinsic alignment model from a lower signal (the lin- 
ear matter power spectrum is lower than the non-linear 
matter power spectrum on small scales). 

4. PHOTOMETRIC REDSHIFT QUALITY 

Photometric redshifts exploit spectral features to use 
images taken with broad wavelength filters to estimate 
redshifts and hence ma ke a three dimensional map of 
the universe. We follow lAmara and Refregierl (|2006[ ) in 
using two parameters to describe photometric redshift 
quality. The first quantifies nearby scatter about the 
true spectroscopic redshift Sz = (7z/{l + z), where cr^ is 
the width of the Gaussian contribution centered on the 
true redshift. The second allows a fraction /cat of catas- 
trophic outliers to exist at a distance ±Az from the cen- 
tral Gaussian distribution. Throu ghout we set = 1 
for simplici ty. See the app endix of lAmara and Refregierl 
(j2006[) and lAbdalla et all (|2007D for more details. 

4.1. Number of redshift bins 

We first consider the effect of the number of red- 
shift bins used for the source galaxies on cosmological 
parameter uncertainties. For this investigation we as- 
sume perfect photometric redshifts by default. It has 
been noted by several a u thors for the s hear s i gnal (GG) 
alone C e.g. IHuI (119 9^; 'Simon et al.' ('2 00l: iMa et all 
(|2006D : lAmara and Rcfregier (2006); Jai n et all (|2007f )) 
that although a small amount of tomographic informa- 
tion greatly improves parameter constraints, continuing 
to subdivide into further redshift bins beyond ~ 3 does 
not significantly reduce uncertainties further. This can 
be seen by the dotted line in Fig. [B) Using about 3 red- 
shift bins brings parameter constraints to within about 
20 per cent of the best obtainable. 

However on including intrinsic alignments we expect 



to need more information to constrain more parameters. 
The amount that the constraints on cosmological param- 
eters are weakened will depend on the number of free pa- 
rameters in our intrinsic alignment modeling. We start 
by considering 5 bins in wavenumber plus an overall am- 
plitude and redshift power law (7 parameters for each 
of II and GI) . This corresponds to the circles marked in 
Fig. [3 

The solid line in Fig.[S]shows that to reduce uncertain- 
ties to 20 per cent of the minimum we now need about 6, 
7 or 5 redshift bins for as, f^m and rig respectively. We 
need 10 redshift bins to bring the FoM to 80 per cent of 
its maximum value. The maximum number of redshift 
bins we have used is 20, but if we were to use more then 
the numbers of redshift bins required would be larger. 
Thus about twice as many (or more) redshift bins are re- 
quired to extract the best possible results from the data 
for the parameters us, flm and twice as many bins are 
required, whereas three times as many are needed for the 
FoM. 

For the above we have assumed perfect photometric 
redshifts. We now consider a photometric redshift error 
of 6z = 0.05. Therefore at the median survey redshift 
the photometric redshift uncertainty is ^2(1 -I- z) 0.1 
so there is potential for effectively about 20 redshift bins 
between < 2; < 2. The non-zero photometric redshift 
error has negligible effect on the lensing alone (GG) curve 
and so we do not plot it. This is in keeping with the fact 
that a small number of redshift bins are sufficient to ex- 
tract most of the information. For the cases in which 
wc have marginalized over intrinsic alignment parame- 
ters then the non-zero photometric redshift error simply 
causes the curves to converge to the best possible FoM 
at a smaller number of redshift bins. This is because the 
best possible FoM is limited by the photometric redshift 
quality, and in effect sets a maximum number of red- 
shift bins, thus increasing the best obtainable parameter 
uncertainty. As an example we show the result for our 
default intrinsic alignment model by the dashed lines in 
Fig. El 

We return to perfect photometric redshifts to check 
the effect of our choice of parametrization of the intrin- 
sic alignments. We considered reducing the number of 
wavenumber bins from 5 to 2 (while keeping just one 
redshift bin). This made negligible difference compared 
to our default parametrization (solid line) so we do not 
show it in the Figure. We next switch to our maximally 
flexible intrinsic alignment model (n = 10, to = 5). The 
dot-dashed lines in Fig. [6] show how much the required 
number of redshift bins is increased. We now need 8,9, 
10, 12 bins to reach 20 per cent of the best uncertainty 
for as, ^m, ng and the FoM respectively, a factor of 3 or 
4 more than that required by lensing alone. 

We repeated all the above results using a different fidu- 
cial model {as = 0.9, = 0.3) and found that these rel- 
ative results changed very little, even though the param- 
eter constraints are tighter for this model with a larger 
signal. We see that the results are similar if using the lin- 
ear matter power spectrum in the linear alignment model 
(light solid line in Fig. [6]). The photometric redshift re- 
quirements are less stringent but this is partly because 
the constraints are already more degraded (see Fig. [5]). 

4.2. Photometric redshift error 
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Fig. 6. — Dependence of cosmological parameter uncertainties on the number of redshift bins into which the source galaxies are 
subdivided. We normalise the performance relative to that for 20 source redshift bins, to show the fraction of the maximum information 
that is obtained. For absolute dark energy Figure of Merit values see the previous figure. We find that intrinsic alignments require more 
redshift bins to achieve the maximum potential. Top left: Uncertainty divided by minimum uncertainty (at 20 lensed galaxy redshift 
bins) for the fiuctuation normalization parameter erg. Top right: Same for the matter density parameter Cm. Bottom left: Same for the 
power spectrum slope parameter ris . Bottom right: Dark energy Figure of Merit, divided by the maximum Figure of Merit (for 20 lensed 
galaxy redshift bins). Dotted lines: Results if intrinsic alignments did not exist. Solid lines: Results for our default intrinsic alignment 
model. Dashed lines: Results using a more realistic photometric redshift error of 5^ = 0.05 (all other lines assume 5^ =0), default intrinsic 
alignment flexibility. Dot-dashed lines: A more flexible intrinsic alignment model. Light solid lines: Using the linear theory matter power 
spectrum in the linear alignment model instead of the non-linear theory matter power spectrum, default intrinsic alignment flexibility. 



Although it is quahtatively helpful to think about the 
number of photometric redshift bins, in practice we will 
have a limited photometric redshift accuracy and will use 
as many redshift bins as computationall y possi ble (or a 
fully three dimensional analysis Heaven s" et al.l ([2006)). 
We therefore repackage our results in terms of the pho- 
tometric redshift accuracy parameter Sz- We also inves- 
tigate the relative importance of the inclusion of the II 
and GI terms, but focus in on the FoM for dark energy. 

The photometric redshift uncertainties expected from 
future surveys depend critically on the number of 
observing filters. Using real data with approxi- 
mately four optical observing ba nds a typica l scat- 
ter az is about 0.05 to 0.1 (jC sabai et al.' '2003'; 
Collister and Laha'j [20041: IPadmanabhan et al. 2005; 



by re r noving physically close pairs (jKing and Schneideil 
20021: iHevmans and Heavens! l2003t iTakada and White! 



Ilbert et al.!l2006i : [Abdalla et al."2007). 

Considering the lensing terms alone, we find again the 
usual result that the constraints are not very sensitive to 
photometric redshift quality (dotted lines in Fig. [7]). 

Ignoring intrinsic-shear correlations for now, we see in 
the left hand panel of Fig. |7| that the II term always 
prefers better photometric redshifts, when the number of 
redshift bins is large i.e. the solid line does not flatten at 
low 5z ■ This makes sense because the correlation length 
for intrinsic alignments (~ 10 Mpc) is small compared to 
the photometric redshift errors. For instance, at z = 1, a 
comoving distance of 10 Mpc corresponds to a 5z = 0.004, 
and a.Sz= 0.05 is 125 Mpc. 

Some works propose to remove the II contamination 



2004! ) ■ instead of marginalising over a model. It is clear 
that for this technique the quality of the photometric 
redshifts is paramount if a large number of pairs are not 
to be rejected. 

To take only a 20 per cent reduction in the FoM 
(relative to the best possible) the lensing alone allow 
6z ^ 0.08, whereas to remove the II term of intrinsic 
alignments requires twice the precision {5z ^ 0.04). The 
left hand panel of Fig. |7| shows that this result is rela- 
tively independent of the model used for intrinsic align- 
ments (solid and dot-dashed lines are similar at low 6z )- 
A very flexible model does very badly for poor photo- 
metric redshift information. As expected, the number 
of redshift bins is most important at high photometric 
redshift accuracy (low Sz). 

Considering just the GI term (and ignoring the II), the 
middle panel of Fig. [7] shows roughly similar results to 
the left hand panel (II but no GI) . The curves are flatter 
at high photometric redshift accuracy (smaller Sz) which 
can be explained by the inclusion of the broad lensing 
kernel in the GI redshift dependence. The number of to- 
mographic redshift bins is not crucial for removing the GI 
term from the GG (the solid, 20 bins, and dashed, 5 bins, 
lines are very similar). However the photometric redshift 
quality is still important. A twenty per cent reduction in 
relative FoM requires approximately 30 per cent smaller 
redshift uncertainties than if no intrinsic alignments are 



11 



GG and II 



GG and Gl 



GG, Gl and II 




GG alone 

- Default 

- 5 source z bins 

■ Maximally flexible 

- Linear model 




Fig. 7. — Relative dark energy Figure of Merit as a function of photometric redshift uncertainty parameter 5z = 0-2/(1 + z). The FoM 
degrades more rapidly with decreasing photometric redshift quality when intrinsic alignments are considered. Left panel: Including just 
the II term with the usual lensing terms i.e. assuming that the Gl contribution is zero, and known to be so. Middle panel: Including just 
the Gl term (with GG) i.e. assuming that the II contribution is zero, and known to be so. Right panel: Both II and Gl terms included. 
Dotted lines: Results using just lensing terms alone (therefore these are the same in every figure), i.e. assuming that the II and Gl terms 
are zero, and known to be so. Solid lines: Using our default intrinsic alignment model parametrization. Dashed lines: Using 5 source 
galaxy photometric redshift bins (the other lines use 20 bins), default intrinsic alignment model. Dot-dashed lines: A more flexible intrinsic 
alignment model (n = 10, m = 5). Light solid lines: Using the linear theory matter power spectrum in the linear alignment model instead 
of the non-linear theory matter power spectrum, default intrinsic alignment flexibility. 

considered. 

Including both the Gl and II terms gives a bigger hit on 
parameter accuracy than only one term alone. This is not 
surprising since we have significantly increased the num- 
ber of free parameters, without adding any more data. 
The results look qualitatively like those for GG-I-II alone, 
but place even more stringent constraints on photometric 
redshift quality. To obtain 80 per cent of the best possi- 
ble FoM, the photometric redshift errors now need to be 
a factor of 3 to 4 better than if considering lensing alone, 
depending on the number of parameters considered in 
the intrinsic alignment analysis. 

Slicing in the opposite direction, if the photometric 
redshift quality is about Sz ~ 0.05, the lensing con- 
straints are degraded by only 10 per cent, relative to that 
for perfect photometric redshifts. Whereas additionally 
considering II, Gl or Gl-f II terms degrades the FoM by 
25, 20 and ~ 40 per cent, relative to that with perfect 
photometric redshifts. 

For all the combinations, we find the results are not 
very sensitive to the number of wavenumber bins used 
in the II and Gl parameterization (n=5 and n=20 bins 
give essentially the same curves, whether m=l or ni=5 
redshift bins for the II and Gl parametrization are used) . 
The light solid line in Fig. [7] shows the result using the 
linear matter power spectrum in the linear alignment 
model. Again we see that the requirements are simi- 
lar but slightly less stringent. We repeated the default 
calculations using an outlier fraction /cat = 0.1 but find 
the results to be negligibly different, i.e. all the FoMs 
are degraded by the same factor on changing Sz- 

We repeated the calculations shown in Fig. [7] but vary- 
ing /cat instead of Sz on the x-axis of the figure. How- 
ever we found that there was little difference between 
the different lines. All dropped to 80 per cent of the 
maximum FoM at /cat ~ 0.1. When intrinsic alignments 
(GG-f GH-II) were included /cat needed to be 10 per cent 
smaller than when intrinsic alignments were assumed to 
be negligible (GG alone). When fewer redshift bins were 
used the requirements on /cat were weaker (5 redshift 
bins instead of 20 allowed the required /cat to increase 
by 30 per cent). This shows that it is useful to carry 
out an analysis with more than 5 redshift bins to make 
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Fig. 8. — The effect on the dark energy figure of merit as priors 
on the intrinsic alignment bin parameter Bij are decreased. Dotted 
line: Results if intrinsic alignments did not exist. Dot-dashed line: 
Results if Gl did not exist. Dashed line: Results if II did not exist. 
Solid line: Results using our fiducial intrinsic alignment model 
including all contributions. 

maximum use of the information available. 

5. EFFECT OF PRIORS 

So far we have placed very wide priors on the pa- 
rameters of the intrinsic alignment model, and assumed 
that the photometric redshift distributions are perfectly 
known. We now consider narrowing the priors on intrin- 
sic alignment parameters and allowing uncertainty in the 
lensed galaxy redshift distributions. Throughout the be- 
low we place the lensed galaxies in 20 photometric red- 
shift bins, assume default photometric redshift parame- 
ters 6z — 0.05, /cat — and default intrinsic alignment 
model with 5 bins in wavenumber, and the same power 
law in redshift for each wavenumber {n = 5, m = 1). 

5.1. Priors on Intrinsic Alignment Model 

We now investigate the impact of tighter constraints on 
the intrinsic alignment mod eling, that could come from 
observations like those of iMandelbaum et al] (120061 ). 
Hira ta et al. (200 7) or advances i n theoretical predic- 
tions (e.g. IHevmans et al.l l2006bl : Ivan den Bosch etldl 
|2002|) . We have previously been using a very wide Gaus- 
sian prior on each bin coefficient Bij of width ABij — 10 
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for numerical stability. Fig. [8] shows how the FoM in- 
creases as this value is decreased. We can see that the 
wide prior we used for the earlier sections was sufficiently 
wide that the self-calibration regime was reached. We 
also see that if the priors on the lA model are sufficiently 
tight then the constraints on cosmology (GG-f Gl-fll) are 
tighter than if intrinsic alignments were negligible (GG). 

As in Fig [3 GGH-Gl gives a lower FoM than GG+II. 
In addition the transition to the self-calibration regime 
occurs at smaller ABij for GG+Gl than for GG-I-II. This 
means that more precise information is required to alle- 
viate the cosmological parameter degradation, compared 
to GG-I-II. The exact transition to self-calibration will 
in general depend on the number of parameters being 
fitted, however for our default model with 7 parameters 
for each of II and GI, priors on the perturbing parame- 
ters Bij need to be less than 10 per cent to help improve 
cosmological parameters constraints. We repeated the 
calculation keeping the overall amplitude and power law 
parameters fixed but the qualitative picture remains the 
same. 

5.2. Priors on redshift distributions 

In practice the redshift distributions of the Icnscd 
galaxies cannot be known perfectly because we only have 
imperfect photometric redshifts. So far we have dealt 
with the accuracy properties of each single photometric 
redshift estimate, in terms of a scatter parameter and 
a catastrophic outlier parameter. For cosmic shear it is 
generally much more important to have accurate knowl- 
edge of these parameters than for these parameters to be 
small. It is difficult to know the values of these parame- 
ters accurately since t hey depend ori uncer tain modeling 
of galaxy formation (jAbdalla et al.l [20071 ) and/or large 
test samples for which spectroscopic redshifts are known. 

Here, we discuss our knowledge of these parame- 
ters by considering our uncertainty on them, and place 
Gaussian priors on their va l ues. This follows previ- 
ous w ork (iHuterer et al.|[2006l:lMa et al.l l2006l: iJain et all 
120071: lAmara and Refregieil I2006D . which has also dis- 
cussed the numbers of spectroscopic redshifts needed to 
obtain these priors. We revert back to wide priors on the 
intrinsic alignment parameters ABij = 10. 

We divide the redshift range < z < 3 up into 30 pho- 
tometric redshift prior bins such that an equal number 
of galaxies fall into each bin. We apply a shift of 2;bias(i) 
to all the spectroscopic redshifts in bin number i. We 
assign the scatter for bin i to be 6z{i){l + z). We then al- 
low these 30 parameters to vary in the Fisher matrix and 
apply priors of width A^^ and Azbias which are the same 
for every photometric redshift prior bin. We then create 
the joint probability distribution P{zs, Zp) and make cuts 
in photometric redshift space, as before. 

This is the same as in lMa et al.l (|2006[ ) except we vary 
Sz instead of cr^ = 5^(1 -I- z) and using binning with equal 
numbers of galaxies in each bin instead of even spacing in 
redshift. The broad results a re similar but m ore accen- 
tuated when we switch to the lMa et aT] (|2006D method. 

The results in Fig. [S] show that in terms of absolute 
FoM, tighter priors are required when intrinsic align- 
ments are considered, if we are to obtain the same FoM 
as in their absence. For example, to obtain an FoM 
of 10 extremely tight priors are required if considering 
GG+GI+II, whereas a wider prior could be used if in- 
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Fig. 9. — The effect of loosening our priors on photometric 
redshift parameters. Previously we effectively assumed zero for 
the uncertainty on photometric the redshift scatter AS^il^ z) and 
on the uncertainty in photometric redshift bias Azbias- Lines as in 
previous figure. 

trinsic alignments did not exist (GG). However, if we are 
considering the relative figure of merit then the results 
are reversed. To obtain 80 per cent of the best possi- 
ble FoM, the presence of intrinsic alignments allows us 
to loosen our requirements on the priors, and thus allow 
smaller spectroscopic redshift training sets. Note that 
this is because the best possible FoM is lower when we 
take into account intrinsic alignments. 

If we had particularly poor prior knowledge of photo- 
metric redshift properties, then the existence of intrin- 
sic alignments actually helps to calibrate the uncertain 
parameters. This is particularly true for the intrinsic- 
intrinsic term. Note that this all assumes our default 
photometric redshift fiducial model in which Sz — 0.05 
and /cat = 0. Therefore if photometric redshifts were 
good, but we did not know it, then intrinsic alignments 
would help. 

For our example with 30 photometric redshift param- 
eter bins the FoM with intrinsic alignments included 
(GG-I-GI-I-II) is better than that if intrinsic alignments 
did not exist for Azbias = ASz > 0.03. However, a 
considerable fraction of the FoM has been lost by this 
sta ge so we hope not to be in this regime. FoUow- 
ing lMa et al.l (|2006f ) we would already have a prior tighter 
than this if we had a number N^p^^ — (tTz/Azbias)^ — 

(('5z(l + 2))/Azbias)^ ^ 10 spectroscopic redshifts in each 
photometric redshift parameter bin. It is likely we will 
already have more spectroscopic redshifts than this, but 
perhaps this feature might become useful for very deep 
surveys, for which spectroscopy will be harder to obtain, 
or as a cross-check. 

6. CONCLUSIONS 

We have compared our fiducial model to current obser- 
vations and an alternative model and found it to be in 
good agreement. We use the linear alignment model as 
our base, and investigate inserting the non-linear matter 
power spectrum, as a way of taki ng into account sma ll 
scale correlations (as suggested in iHirata et al.l ()2007l )). 
We find that this matches better to the HRH* model 
for intrinsic-intrinsic (II) ali gnments, and matches bet - 
ter to the G I detections of iMandelbaum et al.l (120061 ): 
IHirata et al.) ([2007). 

We allow this model freedom as a function of scale and 
redshift since we will not know the intrinsic alignment 
contributions perfectly in time for future dark energy 
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surveys, due to the dependence on baryonic physics. We 
choose as our basehne model one which has five free pa- 
rameters as a function of scale, and the same power law 
function of redshift for each scale. We allow the option 
of adding more redshift variation, and more or less pa- 
rameters as a function of wavenumber. We assume the II 
and GI perturbations about the linear alignment model 
are independent of each other. 

We find that if little freedom is given to the intrinsic 
alignment modeling then the constraints on dark energy 
parameters can actually be very slightly improved, rela- 
tive to those obtained if intrinsic alignments do not exist. 
Essentially the intrinsic alignments are providing another 
cosmological probe. 

The shear-intrinsic correlations are more subtle since 
they contain the lensing efficiency function. It seems 
plausible that cross correlations between a given red- 
shift bin and the various other redshift bins will al- 
low simultaneous constraints on the GI term and cos- 
mology. Thi s is not s o dissi milar to the le nsing ge- 
ometric test fiJain and Ta vlor 20031 : iHu and Jain 200^; 
iBernstein and Jainll2004D T in which the mass of halos is 
jointly constrained with cosmology due to the depen- 
dence of the shear signal with redshift. However, on 
allowing realistic freedom in the intrinsic alignment mod- 
eling and using more realistic photometric redshift uncer- 
tainties the potential is lost, and dark energy constraints 
are degraded by about 50 per cent. 

We show how the number of source galaxy redshift 
bins affects cosmological parameter constraints with and 
without intrinsic alignments. Wc find that on consider- 
ing intrinsic alignments the number of redshift bins must 
be doubled or tripled to obtain the full information avail- 
able. 

We investigate in more detail the effect of photometric 
redshift quality on constraints. We consider how con- 
straints are degraded with respect to the case in which 
perfect photometric redshifts are available. We confirm 
that intrinsic alignments put tighter constraints on pho- 
tometric redshift quality than lensing alone. We find 
roughly equal contributions from each of the II and GI 
terms. The photometric redshift error required is ap- 
proximately three times smaller that using lensing alone, 
to extract the same fraction of the available information. 
It would therefore be most unwise to plan future grav- 
itational lensing experiments to capitalize on the result 
that lensing does not require precise redshift informa - 
tion. The color tomography method of lJain et al.l ()2007l ) 
could likely still be applied if enough observing bands 
were used. 

A significant catastrophic outlier fraction would be bad 
news for methods for removing the intrinsic-intrinsic cor- 
relations by excluding pairs of galaxies from the analysis. 
Statistically it may be known that there is a ten per cent 
chance the true redshift is far from that measured. De- 
ciding whether or not to remove the pair is then more 
unclear, and contamination will leak in. Modeling the 
intrinsic-intrinsic correlations means that statistical in- 
formation on the redshift distribution is sufficient i.e. we 
can use the number as a function of redshift for each bin 
rather than a precise redshift of each galaxy. However 
we do need to use an intrinsic alignment model that can 
e ncapsulate the true i ntrins ic alignment contribution. 

iMandelbaum et all ()2006f ) suggested excluding bright- 



est cluster galaxies (BCGs; which often pass the cuts 
used to select luminous red galaxies, LRGs) from samples 
used to measure cosmic shear, in order to minimize the 
intrinsic terms. However, these galaxies do have smaller 
photometric redshift errors and hence would be useful to 
include in 3-D lensing studies, as well as in the constraint 
of cosmological parameters that we have considered here. 
Perhaps simultaneously fitting for the intrinsic terms as 
well as for the cosmic shear signal as a function of galaxy 
type (plus relative redshift and scale), rather than ex- 
cluding BCGs would be possible. In this paper we have 
ignored the fact that the intrinsic alignment parameters 
depend on galaxy type. 

Obtaining additional information on the intrinsic align- 
ments from other sources would be extremely valuable, 
since the degradation of cosmological parameter con- 
straints is large on using a very flexible model for intrin- 
sic alignments. This could come from better numerical 
simulations or from observations of the galaxy-shear and 
galaxy-ga laxy cross co rrelation functions such as those 
bv MandelLaum et al.l (12006.): Hirata et all ()2007D or as 
suggested bv lHu and JainI ( 20041 ). We find that priors of 
around 10 per cent would be helpful. 

Cosmic shear produces only E-mode (curl-free) distor- 
tions, except on sub-arcminute scales where source clus- 
tering induces a B-mode signal (Schneider et al. 2002). 
The presence of B-mode distortions would indicate that 
residual systematics are present, either from incomplete 
correction of galaxy ellipticities for distortions (associ- 
ated with the atmosphere or telescope) or from intrin- 
sic galaxy alignments. The magnitudes of the B-mode 
power spectra of II and GI are highly model dependent 
- thus, if one can be sure that the data have been ac- 
curately corrected for other systematics, this signal will 
place stringent constaints on which physical models of 
intrinsic alignment are supported by the data. 

To summarise our findings: 

• When constraining the dark energy equation of 
state using a cosmic shear survey, it is essential to 
account for intrinsic-intrinsic galaxy alignments (II 
correlations) and the cross-term between intrinsic 
ellipticity and cosmic shear (GI correlations). Ne- 
glecting these correlations biases the dark energy 
equation of state by ^--^50%. 

• We considered a range of flexibility in the intrin- 
sic alignment modelling and found a wide variation 
in the impact on the degradation in the dark en- 
ergy figure of merit. It will be necessary for future 
observational and theoretical studies to inform this 
choice when making future dark energy constraints. 

• In general when II and GI correlations are ac- 
counted for in cosmological parameter estimation, 
the number of tomographic redshift bins must be a 
factor of at least two more than if only lensing cor- 
relations are present, to obtain 80% of the available 
information. 

• Concerning the dark energy figure of merit, the 
presence of II and GI correlations places a more 
stringent requirement on the accuracy of photomet- 
ric redshifts - a factor of ~ 3 - compared with the 
case where only cosmic shear is present. 
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